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. Abstract 

^ : 

COSY data on the cross section of the reaction pp — > {pp}s^ , where {pp} s is the proton 
pair in the state at small excitation energy E pp = — 3 MeV, recently obtained for beam 
energies 0.5 - 2.0 GeV are analyzed within the one-pion exchange model. The model is based 
on the subprocess ir°p — ► tt°p and final state pp-interaction. A broad maximum observed in the 
+f • energy dependence of the cross section at 0.5 - 1.4 GeV in the forward direction is explained 

. by this model as a dominant contribution of the isospin | in the 7r°p-rescattering. The second 

maximum observed at 2 GeV is underpredicted within the model by one order of magnitude. 
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<n : 

cn 1 Introduction 

<N ; 

Study of the reaction pp — > {pp} s 7r°, where {pp] s is the proton pair (diproton) in the 
1 S state at small excitation energy E pp = — 3 MeV, is motivated by several reasons. 
First, this is the simplest inelastic process in the pp-collision, which can reveal underlying 
dynamics of NN interaction. Second, restriction to only one pp-partial wave (s-wave) 
^ , in the final state considerably simplifies a comparison with theory making it basically 

similar to that for the other simplest reaction of this type, pp — ► dir + . However, while 
for the reaction pp — ► dir + there are a lot of data including spin observables [1], which 
are used to test theoretical models in the GeV region [2, 3], data on the reaction pp — > 
{pp} s TT° above 0.4 GeV were absent until recent measurements at COSY [4, 5]. Third, 
the quasi-binary reaction pp — > {pp} s vr° is very similar kinematicaliy to the reaction 
pp — > dn + , but its dynamics can be essentially different. In fact, quantum numbers of 
the diproton state (J 71 = + , / = 1, S = 0, L = 0) differ from these for the deuteron 
(J 77 = + , I = 0, S = 1,L = 0,2). Therefore, transition matrix elements for these two 
reactions are also different. Using the generalized Pauli principle and angular momentum 
and P-pariry conservation, one can easily find that only negative parity states are allowed 
in the reaction pp — > {pp} s n°. Thus, for the intermediate AY state odd partial waves 
(p-, f-, . . .) are allowed, whereas even waves (s-, d-, . . .) are forbidden. Therefore, at the 
nominal A(1232)-threshold of the reaction YY — > AY, T p = 0.63 GeV, the lowest allowed 
partial wave is the p-wave, which, however, has to be suppressed by the centrifugal barrier. 
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In contrast, in the pp — > dir + reaction both negative and positive parity A — N states are 
allowed. As a consequence, the relative contribution of the A-mechanism to the reaction 
pp — > {pp} s n° is expected to be suppressed as compared to the reaction pp — ■> dn + . This 
argument was applied in Ref. [6] to explain a very small ratio (less of few percents) of 
the spin-singlet to spin-triplet pn-pairs observed in the LAMPF data [7] in the final state 
interaction region of the reaction pp — > pnir + at proton beam energy 0.8 GeV. Obviously, 
this argument is valid for any intermediate N*N- states with other nucleon isobars iV* of 
positive parity. Furthermore, since A— type mechanisms are of long-range type, reduction 
of their contribution would mean that other mechanisms, like iV*-exchanges [8] which 
are more sensitive to short-range NN-dynamics, could be more important in the reaction 
pp — > {pp} s 7T° as compared to the pp — > dn + reaction [9]. 

The cross section of the reaction pp — > {pp} s ir° was measured recently at energy 0.8 
GeV in Ref. [4] and at beam energies 0.5 - 2.0 GeV in Ref. [5]. (For measurements at en- 
ergies below 0.425 GeV see Refs.[10, 11, 12].) At the zero angle, the data [5] show a broad 
maximum in the energy dependence of the cross section at 0.5 -1.4 GeV. This maximum is 
similar in shape and position to the well known A— maximum in the reaction pp — > dn + . 
However, a comparison with the calculation [13] performed within a microscopical model, 
which includes A(1232)-isobar excitation and s-wave 7riV-rescattering, shows very strong 
disagreement between the model and the data obtained at energies 0.5 - 0.9 GeV [5] both 
in the absolute value and shape of energy dependence of the cross section. So, the forward 
cross section measured in Ref. [5] is lower than the calculated one [13] by factor of three 
at 0.6 GeV, where the A— isobar maximum would be expected, whereas at 0.5 and 0.8 
GeV the disagreement is more than one order of magnitude [5]. 
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Figure 1: The OPE mechanism of the reaction pp — > {pp} s 7r°. 

In view of qualitative arguments given above, this disagreement would mean that 
either the model of Ref. [13] is incorrect, or the observed maximum of the cross section 
of the reaction pp — > {pp} s 7r° at 0.5 - 1.4 GeV is of non-A-isobar origin 2 . Here we 
analyse these data employing a simpler model, which includes the subprocess ir°p — > n°p 
and the final state pp( 1 5 , )-interaction (Fig.l). We show that the observed shape of the 

2 It was supposed in Ref. [5] that the observed maximum is caused by the A— isobar contribution, but 
no theoretical calculation were performed to confirm it. 
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peak and, to a large extent, its magnitude are in agreement with the dominance of the 
A(1232)-isobar contribution. 



1.1 The model 

We consider the reaction pp — > {pp} s TT° within the mechanism which corresponds to the 
triangle diagram in Fig. 1. A very similar mechanism was successfully applied for analysis 
of the pp — > dn + reaction in the region of the A(1232)-isobar [14, 15] and at higher energies 
too [14]. 

The amplitude of the reaction pp — > {pp} s 7i° consists of two terms, A = A dir — A exch , 
where A dir is the direct term and A exch is the exchange one. These terms are related to 
one another by permutation of two initial protons. The one-loop integral for the direct 
term A dtr is evaluated very similarly to the OPE-II model of the reaction pd — > {pp} s n 
considered in Ref.[16]. Thus, A dir takes the following form: 

A dir (pi, a u p 2 , a 2 ) = f -^LN pp 2m p F nNN {kl) x (1) 



xS CT 3 CT4 ^(^3^ 4 |00)(l/iia3|ia 1 )^(p,7)^(7r p^7r p), 



here f^NN is the ttNN coupling constant with P^nnI^ — 0.0796, F 7rNN (k 2 7T ) = (A 2 — 
m 2 ^) I (A 2 — k 2 n ) is the nNN form factor, k n is the four-momentum of the virtual 7r-meson, 
irip (m n ) is the nucleon (pion) mass, <7j (i — 1, . . . , 4) is the z-projection of the spin of 
ith proton; A(n°p — > rr°p) is the amplitude of the n°p elastic scattering which is taken 
on-mass-shell; the vector J M is defined by the transition form factors as 



where 



f°° (-)* 
Fo{p,l)= / drrjo(prm (r)exp(-7r), (3) 
Jo 

f°° (-)* 

^io(p,7) = * / dr^iprjipl (r)(l + jr) exp (-jr), (4) 
Jo 

here ji(x) (I = 0, 1) is the spherical Bessel function, ip 1 ^ \r) is the pp-scattering wave 
function that is the solution of the Schrodinger equation at the cms momentum |k| with 
the interaction potential V^So) for the following boundary condition at r — > oo: 

4i-M - ==^- (5) 

Here S is the l So phase shift (for simplicity we omit here the Coulomb interaction, which 
is taken into account in real numerical calculations). In Eq.(l) the combinatorial factor 
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Npp = 2 takes into account identity of two protons. Kinematical variables in Eqs. (2) - 
(4) are defined as 

2 _ ml _ m p Ti _ _ pi 

7 ~ (E 1 /m p ) 2 + E~^Jm~ p " " Pl (^7+^Bi"' P "^7^' U 

where pi and T\ — E\ — m p are the total energy, 3-momentum and kinetic energy of 
the initial proton pi, respectively, in the rest frame of the final diproton. The exchange 
amplitude A exch can be obtained from Eqs.(l)-(6) by interchanging 1^2. 
The OPE cross section of the reaction pp — > {pp} s n° in the cm system is 

da 1 Pf f kmaX ,12 k 1 /• -— 
where A; ma:E is the maximal relative momentum in the final pp-system, related to the 



maximal relative energy E rnax as k max = y E max m p , pi (pf) is the cms momentum in the 

initial (final) state of the pp — > {pp} s ^° reaction, s pp is the squared invariant mass of 
the initial pp-system. The factor | in front of the integral over directions of k takes into 
account identity of two final protons. Keeping only the direct term of Eq.(l), one can 
finally find from Eq. (7) 



d\ I e 24n 2 Pi s 



pp 
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N N pp m p F 7rNN (kl) 



x 



x r°^A;^L={2|J^(p,5)| 2 + |J^(p,5)| 2 }^ (A -A)- (8) 

JO J m 2 + k 2l J dil<f> 

The differential crosss section of the reaction n°p — > 7r°p is taken in Eq. (8) at the squared 
invariant mass of the Tip system, s wp , defined as 

s np = (m n + m p ) 2 + 2T n m p . (9) 

Here T n is the kinetic energy of the final meson ir° in the rest frame of the final diproton. 
If 9 is the angle between the cms momenta of the diproton and the proton pi, which emits 
the virtual pion in the direct OPE diagram in Fig.l, and <j) is the cms scattering angle of 
the 7r°-meson in the process n°(k n ) + p 2 — > p 4 + 7r°(g 7r ), then one can find the following 
relation: 



P2oqo 



+ |P2||q^| cos0 = ^m 2 p + p 2 ^/ml + pj - piPfCosO, (10) 



where the four-momenta of the initial proton p 2 = (P20, P2) and the final 7r°-meson q n 
(go, Qtt) in the cms of the np system can be written as 

P20 = 9 \ (s P n + ™ 2 P - kl), g = „ ~ {s^p + ml- m 2 p ), 



q^rl = V?o -ml, |p 2 | = \jp\o- m l- (11) 
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The squared four-momentum of the intermediate 7r-meson is 



kl = 2m 2 p + PiPf cos 9 - y/ml + piyjM^ + p}, (12) 

where M pp is the mass of the final diproton. One can find from Eqs. (10), (11) and (12) 
that backward n°p scattering (0 = 180°) dominates diproton formation in the forward 
direction (9 = 0°). 

Analysis of the reaction pp — > dir + in the A region, performed in Refs. [15] shows that 
the contribution of the pole diagram with the neutron exchange is small but non-negligible 
and being added to the OPE diagram with ttN rescattering improves the agreement with 
the data. For the reaction pp — > {pp} s n° a similar pole diagram seems to be less important 
and is not taken into account here. The point is that at T p = 0.5 — 2.0 GeV and 9 = 0° 
the pp — > {pp} s vertex in the pole diagram of the reaction pp — > {pp} s TT° involves the high 
momentum component of the wave function ip^\q) at the relative momentum between 
protons q = 0.4 — 0.6 GeV/c, but in contrast to the pn — > d vertex, does not contain 
the D-wave which is important for the pole diagram at large q. Furthermore, the S-wave 
component has a node at q 0.4 GeV/c (see, for example, Ref.[16]). 

2 Numerical results and discussion 

In numerical calculation we used the data on the elementary ttN reactions from SAID [1]. 
The scattering wave function tp^ \r) of the pp system at low energy < 3 MeV is largely 
independent of the NN model and is calculated here using the Reid soft core potential 
plus Coulomb interaction [17]. The calculated forward cross section multiplied by the 
factor 0.45 is shown in Fig. 2. When comparing the dotted and full lines in Fig. 2, one 
can see that the contribution of the exchange term \A exch \ 2 is much less important than 
the direct term \A dir \ 2 at 9 = 0°. Thus, we neglect below the term A exch . As seen from 
Fig. 1, the OPE model is in good agreement with the observed shape of the cross section 
at 0.5 - 1.4 GeV. Note that the form factors F (p, 7) and $i (p, 7) in Eq. (2) are smooth 
functions of the beam energy T p . Therefore, the calculated shape of the pp — > {pp} s ^ 
cross section follows mainly the Independence of the 7r°p — > 7r p cross section at the 
cms angle <fi = 180°. The disagreement in absolute value by factor 0.45 corresponds to a 
typical factor of absorptive distortions in the initial pp state [18]. The distortions are not 
taken into account in the present work in view of their dependence on unknown details 
of the production mechanism, in particular, on off-shell behaviour of the 7r°p-scattering 
amplitude. Furthermore, one should note that when calculating the diagram in Fig.l, we 
factor the amplitude of the elastic 7r°p scattering outside the integral sign. Within this 
approximation, the contribution of intermediate A — N states of positive parity is not 
excluded from this reaction as it should be in an exact OPE amplitude according to the 
discussion in the Introduction. For this reason, this simple model cannot provide a precise 
absolute value of the cross section of the reaction pp — > {pp} s 7T°. The disagreement in the 
absolute value of the cross section can be also related in part to the neglected off-shell 
effects in the n°p — > 7v°p amplitude and contribution of the mesons 77, 7/ and uj. 
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Figure 2: The cms cross sections of the reaction pp — > {pp} s 7i° at 6 = 0° versus beam 
energy. Data are taken from Ref. [12] (A) and [5] (•). The dotted curve presents the 
calculated cross section for the incoherent sum of the direct and exchange terms of the 
OPE model amplitude. Other curves are obtained with the direct term A dir only for the 
isospin term as excluded from (dashed line) and included (full line) in the n°p — > n°p 
amplitude, as explained in the text. All curves are scaled by the factor 0.45. 



In order to exhibit sensitivity of the calculated cross section to the A-isobar contribu- 
tion, one can completely exclude the contribution of the isospin | from the 7T°p scattering. 
The isospin decomposition of the A(ir°p — > ir°p) amplitude is the following: 

A(7r°p->7r°p) = ^(ai+2a3), (13) 

here ai (as) is the amplitude with the total isospin \ (§). The cross section of the 7r°p 
elastic scattering can be written as 

da(n°p -> 7T°p) = i{dcr(7r + p) + da(n~p) - dcr(ir°n -> 7r~p)|, (14) 

where da(n + p), da(ii~p) and da{it Q n — > n~p) are the differential cross section of the n + p 
and iT~p elastic scattering and charge exchange reaction 7i°n — > n~p, respectively. After 
the amplitude a-± is excluded from Eq. (13), the cross section of the 7T°p scattering takes 
the form 

da(n°p — > ix°p) = — ^3da(n^p) — da(ix + p) + 3da(n°n — > 7T~p) j. (15) 
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In order to exclude the term az from the ir°p elastic scattering in calculation of the cross 
section of the reaction pp — > {pp} s 7i° one should substitute Eq. (15) instead of Eq. (14) 
into Eq. (8). When we do so, the absolute value of the calculated cross section of the 
reaction pp — > {pp} s n° at 0.4 - 1.1 GeV diminishes by two orders of magnitude and comes 
in strong contradiction with the data (see dashed line in Since the amplitude aa of the 
n N elastic scattering is dominated by the A(1232)— isobar at v /s^/v ~ 1-15 — 1.35 GeV, 
this analysis shows that the excitation of the A(1232)-isobar dominates in the reaction 
pp — > {pp} s it Q at 0.5 - 1.0 GeV too. On the other hand, since the s— wave intermediate 
A — N state is forbidden, but was not excluded from the reaction amplitude within this 
model, the agreement obtained between the calculated and measured shape of the cross 
section suggests that the absence of this s-state in an exact OPE amplitude would be 
not so crucial for the reaction pp — > {pp} s 7T° at 0.4-1.1 GeV, as might follow from the 
qualitative arguments given in the Introduction. In other words, it would mean that 
the p-wave and higher odd waves are not suppressed drastically by centrifugal barrier 
(perhaps, due to long-range character of the A — N interaction) and make a large enough 
contribution to this reaction. 

Let us make some further comments. Firstly, the second maximum of the forward pp — > 
{pp} s 7T° cross section is observed at 1.97 GeV. The forward pp — > dir + cross section exhibits 
a similar maximum [1]. This peculiarity of the pp — > dn + cross section was interpreted 
in Ref. [14] within the OPE model as a manifestation of heavy nucleon resonances in 
the elastic ttN scattering. One can see from Fig. 2 that the OPE model considerably 
underestimates the magnitude of the observed second maximum in the pp — > {pp} s ^° 
cross section. One may suppose that excitation of heavy A's is not sufficient to explain 
the data at 2 GeV and, therefore, other mechanisms of this reaction like N* exchange 
or Reggeon exchange recently discussed in Ref. [16] make a sizeable contribution in this 
region. To choose between the heavy A— isobars excitation and the N* (or Reggeon) 
exchange mechanism one should measure the ratio of the cross sections pp — > {pp} s ^° 
and pn — > {pp} s ^" [16]. 

Secondly, as can be shown, the present model predicts a smooth increase (5-15%) of 
the differential cross section in forward direction at 9 = 0° — 15°, that is in qualitative 
agreement with the data at 2 GeV, but in disagreement at lower energies. A more detailed 
model, with distortions and explicit A-isobars included, has to be developed to describe 
the angular dependence of the cross section. This kind of model considered in Ref. [13] 
below 0.9 GeV was to some extent successful in this respect, while failed to describe the 
energy dependence. 

3 Conclusion 

Arguments, based on parity and angular momentum conservation, show that the S-wave 
A N-intermediate state is forbidden in the reaction pp — > {pp} s n°, when the final pp-pair 
is produced in the 1 .So-state. The microscopical model [13], which takes into account this 
specific feature of the reaction pp — > {pp\ s n°, since includes explicitly the A-isobar contri- 
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bution via coupled NN— and A^A— channels, is in strong disagreement with the observed 
energy dependence of the recently measured cross section of the reaction pp — > {pp} s n° 
at 0.5 - 1 GeV. On the other hand, a rather simple OPE model developed in the present 
work, which includes the subprocess n°p — > n°p and the final state pp( 1 5 , )-interaction, 
reproduces the observed shape of energy dependence of the cross section of the reaction 
pp — > {pp} s 7T° at 0.5 - 1.4 GeV and to some extent agrees with its absolute value. Thus, 
the OPE model clearly exhibits dominance of the A(1232)-isobar in this region, although 
the angular dependence of the pp — > {pp} s ^° cross section is not described within this 
simple model. One should note, that a quite similar OPE model was recently successfully 
applied to the pd — > {pp} s n reaction in Ref.[16] just in the A-isobar region. Therefore, a 
failure of the model of Ref. [13], most likely, is related not to the A contribution itself, but 
rather caused by interference effects with other terms, for example, with phenomenologi- 
cal heavy meson exchange. More insight into the dynamics of the single pion production 
in pN collision can be gained by further measurement of the reaction pn — > {pp} s ^ at 
the same kinematical conditions. It would be also interesting to get data on the reac- 
tion pp — > {pp} s Tc° at higher excitation energy of the final pp-pair, E pp = 3 — 10 MeV, 
where small components of the pp-wave function start to contribute and allow the S-wave 
intermediate AiV— state. 

I am grateful to V.I. Komarov, A.V. Kulikov and V. Kurbatov for useful comments. 
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